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Abstract. We present a new test of the modified gravity endowed with the Vainshtein 
mechanism with the density profile of a galaxy cluster halo observed through gravitational 
lensing. A scalar degree of freedom in the galileon modified gravity is screened by the Vain- 
shtein mechanism to recover Newtonian gravity in high-density regions, however it might not 
be completely hidden on the outer side of a cluster of galaxies. Then the modified grav- 
ity might yield an observational signature in a surface mass density of a cluster of galaxies 
measured through gravitational lensing, since the scalar field could contribute to the lensing 
potential. We investigate how the transition in the Vainshtein mechanism affects the sur- 
face mass density observed through gravitational lensing, assuming that the density profile 
of a cluster of galaxies follows the original Navarro-Frenk- White (NFW) profile, the gener- 
alized NFW profile and the Einasto profile. We compare the theoretical predictions with 
observational results of the surface mass density reported recently by other researchers. We 
obtain constraints on the amplitude and the typical scale of the transition in the Vainshtein 
mechanism in a subclass of the generalized galileon model. 
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1 Introduction 

Discovery of the late-time accelerated expansion of our universe [1, 2] boosted interests in 
nature of gravity on cosmological scales, which might be a key to explore its origin (e.g., 
Refs. [3-5]). Modification of gravity is an alternative approach to dark energy paradigm 
to explain the cosmic accelerated expansion [3-7]. There have been proposed, for example, 
f(R) model [8-12], Dvali-Gabadadze-Porrati (DGP) model [13], and the galileon model [14]. 
In general, it is very challenging to construct a viable model that explains the late-time 
accelerated expansion of the background universe, satisfying the local gravity constraints 
without theoretical plague simultaneously. 

The DGP model is described in the context of the braneworld scenario, which consists of 
a 3+1-dimensional brane embedded in a 5-dimensional bulk. This model has an interesting 
phenomenology, which yields two branches of the Friedmann equation, a self-accelerating 
branch (sDGP) and a normal branch (nDGP). In the sDGP branch, the expansion of the 
background universe self-accelerates at late times without a cosmological constant, while the 
nDGP branch needs to add a stress-energy component with negative pressure on the brane 
to be consistent with cosmological observations [15, 16]. Unfortunately, the sDGP model 
inevitably gives rise to a ghost [17, 18]. Moreover, the sDGP model is disfavored by the 
cosmological observations [19-21]. 

Inspired by the decoupling limit of the DGP model [17, 22], the galileon gravity theory 
has been studied as a possible alternative to large distance modification of gravity (e.g., [14, 
23-33]). This theory introduces a scalar field invariant under the Galilean shift symmetry 
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du<j> 9u4> + bu in the Minkowski space-time, which keeps the equation of motion being the 
second order differential equation. Although the Lagrangian no longer satisfies the Galilean 
shift symmetry in a curved spacetime, one can construct the generalized galileon model whose 
field equation remains a second order differential equation, which simultaneously admits the 
self-accelerating solution in a FRW universe without a ghost instability [34-36]. 

In general, the gravitational force and the gravitational potential are modified on the 
scales of cosmology as a result of the modification of gravity. However, the solar system 
observations require that Newtonian gravity is recovered on those scales [37]. There have 
been proposed a few successful screening mechanisms depending on modified gravity models, 
e.g., the chameleon mechanism [38-40], the symmetron [41, 42], and the Vainshtein mecha- 
nism [43]. One of the notable features of the galileon gravity is the Vainshtein mechanism, 
which hides the scalar field to recover Newtonian gravity in high-density regions. In the 
Vainshtein mechanism, the self- interaction term like (dcj)) 2 ni(j) plays an important roll. Due 
to the nonlinear interaction, the scalar field is suppressed in high-density regions, hence it is 
screened there. Thanks to this screening mechanism, the galileon gravity models can evade 
the solar-system constraints [14] (c.f. [44]). 

It is worth examining whether the screening mechanism works completely or not around 
high-density regions on cosmological scales. The structure of halos of galaxy and cluster of 
galaxies may be useful to test these modified gravity models. Recently, several researchers 
have investigated such a possibility of testing these modified gravity models on the scales 
of halos of galaxy and cluster of galaxies [45-57]. In [58], the authors have investigated a 
constraint on the Vainshtein mechanism in the context of the massive gravity, using gravita- 
tional lensing and velocity dispersion data from galaxies. In a recent work, we investigated 
how the transition in the Vainshtein mechanism may appear in the sDGP model and the 
galileon model [54]. We found that the circular speed of a test particle in the sDGP and 
the galileon model deviates from that in Newtonian gravity at 10% level on the outer region 
of a halo with the Navarro-Frenk- White (NFW) density profile (e.g., [59, 60]). Thus the 
Vainshtein mechanism might not completely hide the effect of the modified gravity on the 
cluster's scales. 

The purpose of the present paper is to investigate how the transition in the Vainshtein 
mechanism in a subclass of the generalized galileon model appears in an observation of a 
halo density profile measured through the gravitational lensing phenomena. We consider 
a subclass of the most general second-order scalar-tensor theory [34-36, 61], in which the 
scalar field modifies the equation connecting the lensing potential and the matter density. 
In this model, the observational quantities measured through the lensing phenomena are 
contaminated by the scalar field's effect, which might leave testable signatures in observed 
halos. Then, we focus our investigation on the signature of the modified gravity by comparing 
theoretical predictions and observations. 

Many works have been done for testing gravity on cosmological scales using, for example, 
the large scale structures of galaxies and the redshift-space distortions (e.g., [62-67]), and 
the integrated Sachs- Wolfe effect through cross-correlations between the cosmic microwave 
background temperature anisotropies and the galaxy distributions (e.g., [68]). The cluster 
abundance (e.g., [69]) provides a test of gravity on rather smaller scales, but it reflects the 
regime of the linear evolution of density perturbations substantially. On smaller scales, other 
tests of deviation from general relativity have been done [70-74, 78]. A halo of cluster of 
galaxies provides a unique chance to test the gravity theory on the scale between the solar- 
system and the large scale structure of galaxies, where the nonlinear effect of the density 
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perturbations plays an important role. 

In the present paper, we focus our investigation on the surface mass density of cluster 
halo measured with gravitational lensing. Cluster lensing surveys are in progress, which 
is useful to obtain mass distribution of clusters over a wide range of radius by combing 
strong and weak lensing data. Recently, measurements of the surface mass density were 
reported using the gravitational lensing [56, 75-77]. The error of the stacked data is small, 
which is useful to test the modified gravity as will be demonstrated below. Umetsu et al. 
derived a mean surface mass profile of four clusters, A1689, A1703, A370, and C10024+17, 
from Hubble Space Telescope and Subaru images, in the range R = (40 — 2800)/i -1 kpc 
at the mean lensing redshift (z[) = 0.32 [75, 76]. They obtained mass profiles with high- 
precision, by combining independent strong-lensing measurements, weak-lensing distortion, 
and magnification. Oguri et al. presented the mass distribution of a sample of 25 galaxy 
clusters, in the range R = (63 - 5010)/i _1 kpc at the mean lensing redshift (zi) = 0.47 and 
the mean source redshift (z s ) = 1.1 [77]. The cluster sample is based on the Sloan Giant 
Arcs Survey (SGAS) from the Sloan Digital Sky Survey (SDSS). They derived the differential 
surface mass density by combining the strong lensing information from the giant arcs and 
weak lensing measurements from Subaru/Suprime-Cam [79] images. Umetsu et al. and Oguri 
et al. fit the data within the general relativity, and they found that the NFW profile favors 
the data. 

We confront the observed surface mass density of clusters with the theoretical prediction 
of the modified gravity model endowed with the Vainshtein mechanism. In the theoretical 
modeling, we introduce the parameters n and e, respectively, which characterize the am- 
plitude and the typical radius of the transition in the Vainshtein mechanism. We obtain 
constraints on f/, and e. This demonstrates the usefulness of the surface mass density to 
constrain the modified gravity model. 

The structure of this paper is as follows. In section 2, we derive our basic theoretical 
formulas in a subclass of the generalized galileon modified gravity to compare with obser- 
vational quantities from gravitational lensing. In section 3, we scrutinize the theoretical 
behavior of the surface mass density profile depending on familiar density profiles. In section 
4, we present the results of the confrontation between the theoretical predictions and the ob- 
servational results, which derive constraints on the model parameters \i and e characterizing 
the modification of gravity. In section 5, we discuss about the results and our interpreta- 
tions. Section 6 is devoted to summary and conclusions. In appendix A, we summarize the 
definitions of the coefficients in the perturbation equations in section 2. In appendix B, we 
summarize the coefficients in the perturbation equations in the original galileon model. 

Throughout this paper, we use units in which the speed of light equals unity, and 
we follow the convention (—,+,+,+). We use the reduced Planck mass Mpi, defined by 
Mpi = 1 / V 8irG with Newton's gravitational constant G. We adopt the Hubble constant 
Hq = 100 /ikm/s/Mpc with h = 0.702, and the matter density parameter at present = 
0.275 [80]. 

2 Formulation 

2.1 Modified gravity model 

The most general second-order scalar-tensor theory was derived by Horndeski [34] for the 
first time, which was recently rediscovered by Deffayet et al. [35] as the most generalized 
galileon theory, which contains four arbitrary functions of the functions of (j> and X = 
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— fJi (pV l ,(j)/2. The Vainshtein mechanism in the generalized galileon theory is discussed 
recently [44], which clarified variety of the solutions of the scalar field in the generalized 
theory. We here consider a subclass of the general second-order scalar-tensor theories in a 
curved spacetime, which is nonminimally coupled to gravity with the action [61], 



S = I d A x^—g 



^F(<P)R + K{<f>, X) - G(<f>, X)D<f> + £ r 



(2.1) 



where K((j), X) and G(4>, X) are the arbitrary functions of <f> and X, F((f>) is the function 
of <j), and C m is the matter Lagrangian. We assume that the matter fields do not have 
direct couplings with the field <j). This corresponds to the Lagrangian with G4 = F ((/))/ 2 
and G5 = in the most general galileon model [36]. The action reduces to the kinetic 
gravity braiding model, for the choice F(<p) = M™ [29, 33]. The original galileon model is 
reproduced by choosing F(<p) = Mp V K{4>,X) = —X and G((j>,X) = (r 2 /M P1 )X, where r c 
is the parameter [29]. 

In this subsection, we summarize the perturbation equations for gravity and the scalar 
field in the cosmological background. In reference [44], the Vainshtein mechanism in the most 
general second-order scalar-tensor theory was investigated. The model in the present paper 
(2.1), a subclass of the most general second-order scalar-tensor theory, was also investigated 
therein. Following [44], we choose the Newtonian gauge, 

ds 2 = -(1 + 2V(t,x))dt 2 + a(t) 2 (l + 2$(t,x))5 ij dx i dx j , (2.2) 

where a(t) is the scale factor. Within the subhorizon scales with the quasi-static approxima- 
tion, we have the following perturbed equations, 

*±$ = -4 !K G8p + A<p, (2.3) 

$ + ^ = -av?, (2.4) 

and 

^'(^ = -*■«*>. < 2 - 5 > 

where y(x) denotes the perturbation of scalar field defined by </>(i, x) = 4>(t)(l + y(x)), 
A represents the Laplace differentiation operator, 5p is the perturbed matter density, the 
Newton's gravitational constant is defined by G = l/(8nF (((>)), and the coefficients a, £, £, 
and A 2 depend on the functions F(cfi), K(cp,X), and G((p,X), whose explicit expressions are 
summarized in appendix A. 

The above equations are derived as follows (see [44] for details). We consider the grav- 
itational and scalar fields on subhorizon scales sourced by a non-relativistic matter density 
perturbation bp. Then, we may ignore time derivative in the field equations, while keeping 
spatial derivatives. Denoting the quantities for perturbations \E', $ and <p by V, we keep all 
the terms (d 2 V) n with n > 1, because the L 2 (t)d 2 V could be large on small scales, where 
L(t) is a typical length scale associated with the background evolution which may be as large 
as the Hubble radius. Within this framework, it has been shown that the terms (d 2 V) n with 
n > 2 do not appear in the gravitational field equations (2.3) and (2.4) for our model (2.1). 
But only the terms of (d 2 V) 2 appear in the scalar field equation (2.5). 
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Note that the coefficients a, £, £, and A 2 are determined by the background field evolu- 
tion depending on F ((/>), K((j>,X), and G(4>,X). The coefficients may take any value because 
F((p), K((fi,X), and G(<fi,X) are arbitrary functions of 4> and X. We should choose F(ip), 
K(4>,X), and G(4>,X) requiring the accelerated expansion of the universe and satisfying the 
conditions of no ghost nor instability which will limit the arbitrary functions. However, 
we here assume that a, £, £, and A 2 may take any value reflecting the arbitrariness of the 
functions in the model, and consider the constraint from observation through gravitational 
lensing. 

Combining (2.3) and (2.4), we have 

, ( *-$ \ a + 2£A 



In the spherically symmetric case, (2.5) yields 



dip _ r f / 8G\ 2 (M{ 



r 



(2.7) 



where we defined the enclosed mass M(r) = 4-7T J Q r dr'r' 2 5p(r'). In the spherically symmetric 
case, (2.7) leads to 



^ _ fog) _ » ! _ J 1+ _ ^«*tS. (2.8) 

a 2 r 2 dr \ dr y 4A 2 ^ V H y r 3 v/l + 8GA2 f M ( r ) 

As we will show in the below (see also [44]), the lensing signature under the influence 
of the Vainshtein mechanism in our model (2.1) is described by the two combination of the 
parameters A 2 £ and (a + 2£)£/2, on which we focus our investigation. In the original galileon 
model, we have non-zero values of A 2 £ and (a + 2£)£/2 (see appendix B). On the other hand, 
in the model of the sDGP model, we have a = —1 and £ = 1/2, then (a + 2£)£/2 reduces to 
zero, while A 2 £ takes a non-zero value. 1 Therefore, in the sDGP model, the relation between 
the lensing potential (<£ — ^)/2 and the matter density perturbation (2.6) is the same as that 
in the general relativity. In this case, our method cannot put a constraint because there is 
no effect of the modified gravity on the gravitational lensing in such a class of models. 

In order to obtain the closed system of equations, we need another equation for the 
matter density perturbation 5p. Because we are considering the model in which the matter is 
minimally coupled with the scalar field, the matter component follows the usual equation un- 
der the influence of the gravitational potential . Instead of specifying the coupled equations 
of bp and we simply assume the NFW density profile or other well-known profiles for the 
matter density perturbation bp (cf. [54]). From N-body simulations in reference [46], it is sug- 
gested that the NFW profile describes the density profile of halos in the DGP model though 
the parameters of the profiles are scaled. This is a supporting evidence of our assumption in 
the present paper. 

2.2 Gravitational lensing in modified gravity 

Now let us consider observational quantities observed through gravitational lensing. The light 
propagation in the Newtonian gauge (2.2) is described by the lensing potential ($ — \£)/2. 



1 In the spatially flat sDGP model, we have A 2 £ = 2(r c /3/3) 2 , where r c — 1/[(1 — Q m )Ho] and /3 
l-2Hr c (l + H/3H 2 ). 
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For example, the convergence k(x) °f the gravitational lensing is given by (e.g., [81, 82]) 

1 rx 

K — 

2 



d , Ik{x ~ X')fK(x') A (2D) ( q 

o x Mx) [ 



(2.9) 



where fxix) is the comoving angular diameter distance and A^ 2D ^ is the comoving two 
dimensional Laplacian. Note that Jk{x) = X f° r the spatially flat universe. The shear can 
be written in a similar form with (2.9). It is worthy to note that any quantity of gravitational 
lensing is described by the lensing potential. 

Using the thin lens approximation, (2.9) and (2.6) yield 



fi<(xs ~ Xl)/k(xl) 



Ys 



4irGp(r') 



a + 2£ A 



at, 



(2.10) 



where xs is the comoving distance between the observer and the source, Xh is the comoving 
radial distance between the observer and the lens object, and <2l is the scale factor at which 
we observe the lens object. We define the observed surface mass density, which is obtained 
through the gravitational lensing phenomena, 



Xs 



d(aLX f ) 



p{r>) 



a + 2£ A 
8vrG 'aJ iP 



where S cr j t is defined by 



1 



J crit 



(2.11) 



(2.12) 



4vrG /x(xl)/x(xs - Xl)«l ' 

We introduce the spatial coordinate being the physical coordinate, whose origin is lo- 
cated at the center of the lens object by Z = ol(x — Xl) and r± = a^fKiXh)^-, where 6 is 
the polar angle of the polar axis connecting the observer and the lens object. Then, using 
(2.8), the observed surface mass density (2.11) is written as 



£s(r±) 



dZ 



p{r) 



3(a + 2Q 
32vrGA 2 



1 



8GA 2 CM(r) 



1 + 



(a + 2£)C 4vrr 3 p(r) - 3M(r) 



8tt 



-i . 8G\ 2 (M(r) 
1 H ^3 



(2.13) 



where r 



+ Z*. 



2.3 Parametrization of the modified gravity 

To investigate the constraint on the modified gravity model on the halo scales, we introduce 
the parameters p and e as follows, which characterize the modification of gravity, 

(a + 2£)C 



/' 



#o 2 A 2 C 



2 V °' 

With these two parameters the surface mass density (2.13) is written as 



(2.14) 



E s (r ± ) = 2 



dZ 



p(r) 



16vrGe 2 



1 



1 + 



8Ge 2 M(r) \ p 4irr 3 p(r) - 3M(r) 



# 2 r 3 



+ 



47T r 3 



1 + 



8Ge 2 M(r) 



(2.15) 
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The physical meaning of fi and e is understood as follows. In the limit of large r, (2.8) and 
(2.6) reduce to 



A 3G(M(r) 
—ip ~ 



and 



A 



2 V 2 



4vrG cfr 5p, 



(2.16) 



(2.17) 



respectively, where G c r = G(l + fx) is the effective gravitational constant in the linearized 
limit. Thus is the amplitude of the modification of gravity in the outer region of a halo in 
the linearized regime. On the other hand, in the limit of small r, (2.7) gives 



dr 



4A 2 



1 l8G\ 2 (M(r) GMir) 

<C 7, ~ 



d^ 

dr 



d&_ 

dr 



(2.18) 



This means that the effect of the scalar field is screened due to the nonlinear interaction of 
the scalar field in the limit of small r, which is the Vainshtein mechanism. Hence, the basic 
equations (2.3), (2.4) and (2.6) reduce to those of general relativity when the scalar field is 
neglected. 

Thus, the modified gravity effect appears in the linear regime at large radii, while the 
modified gravity effect is screened due to the nonlinear effect at small radii. The transition 
radius between the two regime is the Vainshtein radius ry, which we defined by (cf. [33, 54]) 



ry = [8G\ 2 (M V 



I V3 



8Ge 2 M v 



1/3 



(2.19) 



where M v i r is the virial mass of a cluster (see below). Therefore, e is the parameter of the 
Vainshtein radius ry. In summary, Newtonian gravity is recovered due to the Vainshtein 
mechanism for r <C ry, while the gravity is modified for r ^ ry. The Vainshtein radius is 
rewritten as 



ry 



13.4e 2/3 



1O 15 M 



1/3 



/T^Mpc. 



(2.20) 



In the limit \i — > or e — > oo, Newtonian gravity is reproduced on all scales. The original 
galileon model corresponds to [i = 0.26 and e = 0.53 (/i = 0.19 and e = 0.43) at the redshift 
0.32 (0.47), respectively, which is the mean redshift of clusters of the observational data used 
below. 

Figures 1 and 2 demonstrate typical behavior of the surface mass density (left panel) 
and the logarithmic slope of the surface mass density (right panel) as a function of r±. The 
curves are the theoretical model parametrized by /x and e, while the points with the error 
bars show the observational data by Umetsu et al. in [75, 76], which is obtained by averaging 
over four massive clusters. In figure 1, the theoretical curves assume the different values of \i 
with fixing e and the other parameters for the halo density profile, as described in the panels, 
where we assumed the NFW density profile (see the next section). This figure demonstrates 
how the theoretical curve depends on /i. As r± becomes large, the amplitude of Ss(r^) is 
enhanced for positive fx, while it is suppressed for negative [i. The (black) solid curve is 
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Newtonian gravity, while the (magenta) dot and long-dashed curve is the original galileon 
model. Figure 2 is the same as figure 1 but with adopting the different theoretical values. 
Here the different values of e are assumed while \x and the other parameters for the density 
profile are fixed, as described in this figure. As e describes the Vainshtein radius ry, within 
which Newtonian gravity is recovered, Sg(rj_) deviates from Newtonian gravity even at small 
radii for small e (small ry), while £s( r ±) deviates from Newtonian gravity only at large radii 
for large e (large ry). For the case \x < 0, £s( r J_) deviates from the line of \i = in the 
opposite side of the case \x > 0. 



T 




log r ± (h 'kpc) log r_(h ! kpc) 

Figure 1. Surface mass density Sg(r^) (left panel) and the logarithmic slope d\nT,s/dlnr± (right 
panel) as function of r±. The data with the error bar is from Umetsu et al. [75, 76], while the curves 
are the theoretical modeling. The (magenta) dot and long-dash curve is the original galileon model 
(/i = 0.26, e = 0.53), but the other curves assume the same value of e — 0.1 but the different values 
of fj, = 1 (green dashed curve), (black solid curve), and —0.5 (blue dot and short-dashed curve), 
respectively. Note that fi = is Newtonian gravity. We here adopted the NFW profile with fixing 
M vir = 1.5 x 1O 15 M //i and c vir = 7.7. 



3 Halo density profiles and differential surface mass density 

The halo density profile has been investigated as an important consequence of the cold dark 
matter paradigm, which is a key component of the universe. Lots of works have been done 
on the halo density profile with cosmological N-body simulations (e.g. [83-91]). The NFW 
universal profile [59, 60] is reported to work well, but the generalized NFW (gNFW) profile 
has also been studied [92-94]. Some studies indicate that the Einasto profile [95] better fits 
the inner cusps [96, 97, 99]. 

For predicting the surface mass density £s(r_i_), we nee d to determine the halo density 
profile p(r). In the present paper, we adopt the halo density profiles described in this section, 2 
whose validity is suggested from N-body simulations. The halo density profiles could be 
affected by the modification of gravity, however, we adopt the same profile irrespective of 
gravity model, whose validity is partially suggested by N-body simulations for the DGP 
model and f{R) model [46, 51]. In this section, after briefly reviewing cold dark matter 

2 We only consider the contribution of the cluster halo profile itself (1-halo term) to Es(r^), and we neglect 
the neighboring halos (2-halo term), which is only important at large r± [94]. 
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Figure 2. Same figure as figure 1 but with different theoretical models. The black solid curve is 
Newtonian gravity fi = 0, the other curves assume the same value of fi = 1 .0 and the different values 
of e = O.f (green dashed curve), 0.5 (blue dot short-dashed), and 1.0 (magenta dot long-dashed), 
respectively. Note that the curves approach Newtonian gravity as e becomes large. We adopted the 
same NFW profile as figure 1, M vir = 1.5 x 10 15 M Q //i and c vir = 7.7. 




Figure 3. Surface mass density £s(?"±) (left panel) and the logarithmic slope d\nT,s/dlnr± (right 
panel) as function of r±. The data with the error bar is from Umetsu et al. [75, 76], while the curves 
are the theoretical models assuming Newtonian gravity ^ = and the NFW profile with c V j r = 7.7 
and M vir = 1.0 x 1O 16 M (green dashed), 1.5 x W 15 M Q /h (black solid), and 1.0 x 1O 15 M (blue 
dot-dashed), respectively. 



(CDM) halo density profiles (see e.g., [97]), we demonstrate how the observed surface mass 
density depends on the modified gravity. This indicates the potential ability of testing the 
gravity theory with the halo density profile. 

3.1 NFW profile 

We first consider the NFW profile [59, 60] 

[r/r s )(l + r/r s ) z 
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2 3 2 3 

log r 1 (h _1 kpc) log r_(h _1 kpc) 



Figure 4. Same figure as figure 3 but with different values c V i r = 20 (green dashed curve), c V i r = 7.7 
(black solid curve), and c v i r = 1.0 (blue dot-dashed curve), with M vir = 1.5 x 10 15 M©//i fixed. Here 
we assumed Newtonian gravity p = and the NFW profile. 




log r ± (h L kpe) log r_(h ^pc) 

Figure 5. Same figure as figure 3 but for the gNFW profile with 7 S = 1.2 (green dashed curve), 
7,, = 1 (black solid curve), and 7 S = 0.7 (blue dot-dashed curve), respectively, with 7; = 3 fixed. Note 
that 7,5 = 1 is the NFW profile. The other parameters are fixed as M vir = 1.5 x 1O 15 M //i, c v i r = 7.7, 
and p = 0. 



where p s = 4p(r s ) is the characteristic density, r s is the characteristic radius where the 
slope of the density profile changes. We define the virial cluster mass and the concentration 
parameter by 

M vir = ^MsQA™ (3 . 2) 

Cvir = j (3-3) 

r s 

respectively, which can be used as the parameters of the NFW profile instead of p s and r s , 
where r v ; r is the virial radius, A vrr is the virial overdensity, and p C i(zi) = p CT fiH 2 {zi) / Hq is 
the critical density. We take the value A v j r = 120. Note that the virial overdensity and the 
critical density depend on the cosmological background, which we do not specify explicitly. 
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Figure 6. Same figure as figure 3 but for the gNFW profile with 7/ = 3.5 (green dashed curve), 
7; = 3 (black solid curve), and ji = 2.5 (blue dot-dashed curve), respectively, with j s = 1 fixed. Note 
that 7; = 3 is the NFW profile. The other parameters are fixed as M v i r = 1.5 x lO 15 M /ft,, c v h- = 7.7, 
and [i = 0. 




Figure 7. Same figure as figure 3 but for the Einasto profile with T = 0.3 (green dashed curve), 
T = 0.23 (blue dot short-dashed curve), and T = 0.17 (magenta dot long-dashed curve), and with 
r_2 = 0.26/i _1 Mpc fixed. The mass parameter is fixed M vir = 1.5 x 10 15 M Q //i, and we assumed 
Newtonian gravity (x = 0. For comparison, we also plot the NFW profile (black solid curve). 



However, when M v ; r and c v i r are marginalized over, the dependence on the cosmological 
background is not necessarily specified because it is absorbed by the redefinition of M v ; r and 

Figure 3 shows the surface mass density Xs( r ±) (left panel) and its logarithmic slope 
<iln T,s/d In r± (right panel), where we adopted the different values of M v { r with the other 
parameters fixed, which are described therein. Here the concentration parameter is fixed 
Cvir = 7.7. Note that /U = means Newtonian gravity. M v j r changes the amplitude of 
Es(n_)- Figure 4 is the same figure as 3, but with adopting the different values of c v \ r and 
fixing the other parameter. Here we fixed the virial mass as M vir = 1.5 x 10 15 M o //i. The 
slope of the surface mass density significantly depends on c v i r . 
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Figure 8. Same figure as figure 3 but for the Einasto profile wit r_2 = 0.1ft. _1 Mpc (green dashed 
curve), 0.22/i _1 Mpc (blue dot short-dashed curve), and l.O/i Mpc (magenta dot long-dashed curve), 
respectively, with V = 0.23 fixed. The mass parameter is fixed M v ; r = 1.5 x 10 15 M Q /ft,, and we 
assumed Newtonian gravity fx = 0. The black solid curve is the NFW profile. 



As we will describe in the next section, we compare the theoretical model with obser- 
vations by introducing the chi square statistics (see expression (4.1)). Within Newtonian 
gravity, the best-fit value of the chi-squared is 5.8 for 13 degrees of freedom (DOF), which 
we find for M Y u = 1.6 x 10 15 M Q //i and c v i r = 7.9. (see the top line of table 1). 

3.2 gNFW profile 

The generalized parametrization of the NFW model can be written in the form (e.g., [92-94]) 

P(r) = (r/r s ps(l+r/r s )^- (3 ' 4) 
The logarithmic slope of this density profile is defined 

_ d\np(r) 

l3D(r) = —^r, (3.5) 

which reduces to 7 S for r r s , and to 7/ for r 3> r s . The NFW profile is reproduced 
for (7 S , 7/) = (1,3). Cosmological N-body simulations indicate that an inner logarithmic 
slope of the density profile 73i)(r) ~ 1-2 and an asymptotic outer slope 73d(^) ~ 2.5 [97]. We 
introduce the radius r_2 at which the outer slope is isothermal, i.e., 73r>(r_2) = 2. For the 
gNFW profile, r_2 = r s (2 — ^ s )/{li — 2) and the corresponding concentration parameter is 
c-2 = r vir /r_ 2 = c vir (7/ - 2)/(2 - 7 S ). 

Figure 5 shows Sg(rj_) and d In Ss /d In r± adopting different values of 7 S with the other 
parameters fixed as 7; = 3.0, M v j r = 1.5 x 10 15 M Q //i, c v ; r = 7.7, and /i = 0. The effect of 
changing 7 S appears only inside halo. Figure 6 is the same as figure 5, but with adopting the 
different values of 72 with the other parameters fixed as 7 S = 1.0, M vir = 1.5 x 10 15 M o //i, 
Cvir = 7.7, and /i = 0. The behavior of Ss( r _i_) is similar to that of figure 4, which causes a 
degeneracy in the parameter space. We find the best-fit value of the chi-squared x 2 /DOF = 
5.5/13 within Newtonian gravity for (7.5,7;) = (0.7,2.8), using the surface mass density data 
in [76]. This value of the chi-squared is slightly better than the best-fit value x 2 /DOF = 
5.8/13 for the NFW profile with (7 S ,7/) = (1,3). The best-fit value for the gNFW profile 
with (7*;, 7;) = (1.1,3.1) is x 2 /DOF = 6.2/13 within Newtonian gravity (see table 1). 
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Figure 9. Differential surface mass density A£ + (rj_) in as a function of r±. The data with the 
error bar is the result by Oguri et al. [77]. In the panel (a), the curves adopt the different values 
of fx = 1.0 (green dashed curve), fx = —0.5 (blue dot short-dashed curve), and fx = (black solid 
curve), respectively, with e = 0.1 fixed. In this panel the NFW profile with M vir = 4.6 x 10 14 M Q //i 
and c v ir = 5.8 is adopted. The panel (b) is the same as the panel (a), but we adopted the different 
values of e = 0.1 (green dashed curve), e = 0.5 (blue dot short-dashed curve), and e = 1.0 (magenta 
long-dashed curve), respectively, with fx — 1.0 fixed. The other parameters of M v j r and c v i r are the 
same as those of the panel (a). The panel (c) is the same as the panel (a), but we adopted the 
different values of M v i r = 1.0 x 10 16 M Q (green dashed curve), 4.6 x 10 14 M Q //i (black solid curve), 
and 1.0 x 10 14 M Q (blue dot-dashed curve), respectively, for the NFW profile. The other parameters 
are fixed as fx = and c V i r = 5.8. The panel (d) is the same as the panel (c), but adopted the different 
values of c V i r = 20 (green dashed curve), 5.8 (black solid curve), and 1.0 (blue dot-dashed curve), 
respectively. The other parameters are fixed as M v - U = 4.6 x 1O 14 M0 and fx = 0.0. 

3.3 Einasto profile 

Finally, we consider the Einasto profile [95] 



p(r) = p_ 2 exp I -- 



r 
r-2 



(3.6) 
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Figure 10. Same figure as figure 9 but with different theoretical models. The panels (a) and (b) 
assume the gNFW profile, and (c) and (d) do the Einasto profile, respectively. The panel (a) adopted 
the gNFW profile with 7 S = 1.2 (green dashed curve), j s = 0.7 (blue dot-dashed curve), and % = 1 
(black solid curve), respectively. The other parameters are fixed as 7; = 3, M v - n = 4.6 x 1O 14 M //i, 
Cyir = 5.8, and fi = 0. In this panel 7 S = 1 is equivalent to the NFW profile. The panel (b) is the same 
as the panel (a), but with 7; = 3.5 (green dashed curve), 7; = 2.5 (blue dot-dashed curve), and 7; = 3, 
respectively. In this panel, we fixed j s = 1, then 7; = 3 is equivalent to the NFW profile. The panel 
(c) adopted the Einasto profile, with V = 0.3 (green dashed curve), T = 0.23 (blue dot short-dashed 
curve), and V = 0.17 (magenta dot long-dashed curve), respectively. The other parameters are fixed 
r_ 2 = 0.22/i- 1 Mpc, M vir = 4.6 x 10 14 M Q //i, and fx = 0. The panel (d) is the same as the panel (c), 
but with r_ 2 = 0.1/i _1 Mpc (green dashed curve), 0.22ft -1 Mpc (blue dot short-dashed curve), and 
1.0/i _1 Mpc (magenta dot long-dashed curve), respectively, while T = 0.23 fixed. 

where r_2 an d p~2 are the radius and the density at which p(r) oc r~ 2 , which means 
73D(r_2) = 2 and p-2 = /j(r_2). The authors in [88] claimed that CDM halos can be 
more properly described by the Einasto profile than the NFW profile, using Millennium sim- 
ulation. They also claimed that the best-fit value of T increases gradually with the increase 
of the virial mass, from T ~ 0.16 for galaxy halos to T ~ 0.3 for the most massive clusters. 
Figure 7 shows Sg(rj_) and dlnSg/dlnrj_ with adopting different values of T but with 
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the other parameters fixed r_2 = 0.26/i~ 1 Mpc, M V - 1T = 1.5 x 10 15 M Q /h, and \i = 0. T changes 
the slope and the amplitude of Sg(r^) at large radii. Figure 8 is the same as figure 7 but 
with varying r_2 with the other parameters fixed T = 0.23, M v i r = 1.5 x 10 15 Mq//i, and 
fj, = 0. The slope and the amplitude of Xs( r _i_) strongly depend on r_2- 

We obtained the best-fit value of the chi-squared 9.9, 7.2, and 9.6 for 13 DOF for the 
Einasto profile with T = 0.17, 0.23, and 0.3, respectively (see table 1). Thus the Einasto 
profile does not fit the observational data better than the NFW profile or the gNFW profile 
within Newtonian gravity. 

3.4 Differential surface mass density 

Oguri et al. [77] obtained the differential surface mass density, which is defined by 

A£+(r ± ) = £ cri t5> + (r±), (3.7) 
where £ cr i t is defined by (2.12), g+(r±) is the reduced shear 

9+ (r±) = T3 ^ iy , (3.8) 

and the tangential shear is defined by 

/ \ — / , \ / \ S s (< r±) - S s (r ± ) 
7+(r_L) = «(< r±) - K{r±) = , (3.9) 

^crit 

with 

Ss(< r ± ) = \ dr' ± r' ± X s (r' ± ). (3.10) 
r ± Jo 

In figures 9 and 10, we demonstrate the behavior of the differential surface mass density 
AS + (ri) as a function of rj_, comparing with the data in [77]. The panels (a), (b), (c) and 
(d) of figure 9 show how AS + (rj_) depends on fi, e, M v j r and c v i r , respectively, adopting the 
NFW profile. The panels (a) and (b) of figure 10 show the dependence on 7 S and 7; in the 
gNFW profile, while the panels (c) and (d) of figure 10 show the dependence on T and r_2 in 
the Einasto profile, respectively. The behaviors are very similar to those of the surface mass 
density described in the above subsections. 

Similar to the surface mass density, we compare the theoretical prediction with the 
observational data by introducing the chi-squared, (4.8). We obtained the best-fit values of 
chi-squared for the NFW profile, the gNFW profile and the Einasto profile within Newtonian 
gravity as well as allowing the modified gravity, which are summarized in table 2. We find the 
best-fit value of the chi-squared x 2 /DOF = 12.9/17 for the NFW profile within Newtonian 
gravity, with the best-fit parameter M v - ir = 6.6 x 1O 14 M0//i and c V i r = 6.1. This is consistent 
with the result by Oguri et al. [77]. The best-fit parameter is given M v ir = 4.6 x 10 14 M Q //i 
and c v i r = 5.7, when we adopt the same definition for the virial mass M v ; r = 4TTr^ il A(z)p rn (z). 

The result is similar to that using the surface mass density in the point that the Einasto 
profile with large value of T does not better fit the observational data than the NFW profile 
and the gNFW profile within Newtonian gravity. Thus the Einasto profile with T = 0.3 is 
not favored within Newtonian gravity. 
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model 


Newtonian 


gravity 




Modified 


gravity 




{is, 71 ) 


M vir [M Q ] 


Cvir(C-2) 


XrR 


M vir [M ] 


Cvir(C-2) 


/ v 1V1 Kjr 


NFW 


(1, 3) 


2.2 x 10 15 


7.9 


5.8 


2.2 x 10 15 


8.0 


5.7 


gNFW 


(0.7, 2.8) 


2.2 x 10 15 


13.9(8.6) 


5.5 


2.2 x 10 15 


13.9(8.6) 


5.5 


gNFW 


(1.1, 3.1) 


2.2 x 10 15 


6.2(7.6) 


6.2 


2.2 x 10 15 


6.3(7.7) 


6.0 


r 


A/ vir [M ] 


r-2 


XgR 


M vir [M ] 


r-2 


Xmg 


Einasto 


0.17 


2.5 x 10 15 


0.30 


9.9 


2.9 x 10 15 


0.35 


8.7 


Einasto 


0.23 


2.2 x 10 15 


0.28 


7.2 


2.1 x 10 15 


0.26 


6.5 


Einasto 


0.3 


1.9 x 10 15 


0.26 


9.6 


1.6 x 10 15 


0.22 


5.6 



Table 1. Best-fit value of the chi-squared of the surface mass density £s(?"-l) for the NFW profile, the 
gNFW profile and the Einasto profile. Here we fixed the values of 7 S and ji for the gNFW profile, and 
r for the Einasto profile, as described in the table. Xgr ^ s the best-fit value for each halo model when 
Newtonian gravity is assumed, and M v - lr and c vu - or r_2 (in unit of /i~ 1 Mpc) yield the best-fit value. 
Xmg * s t ne best-fit value allowing the modified gravity model. The number of degrees of freedom 
(DOF) is 13 and 11 for Xgr and Xmc respectively. 
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Newtonian 


gravity 




Modified 


gravity 




(7s, 7/) 


M vir [M ] 


Cvir(C-2) 


Xgr 


M vir [M ] 


Cvir(c_2) 


Xmg 


NFW 


(1, 3) 


6.8 x 10 14 


6.1 


12.9 


5.9 x 10 15 


6.6 


11.7 


gNFW 


(0.7, 2.8) 


6.8 x 10 14 


10.9(6.7) 


13.1 


6.1 x 10 15 


11.7(7.2) 


12.2 


gNFW 


(1.1, 3.1) 


6.8 x 10 14 


4.8(5.9) 


12.8 


5.9 x 10 15 


5.3(6.5) 


11.5 


r 


A/ vir [M ] 


r-2 


Xgr 


M vir [M ] 


r-2 


Xmg 


Einasto 


0.17 


7.0 x 10 14 


0.23 


12.4 


6.6 x 10 15 


0.21 


12.2 


Einasto 


0.23 


6.8 x 10 14 


0.23 


14.0 


5.4 x 10 15 


0.18 


11.4 


Einasto 


0.3 


6.7 x 10 14 


0.23 


19.7 


4.7 x 10 15 


0.17 


11.8 



Table 2. Same table as table 1 but for the differential surface mass density A£ + (rj_). The number 
of DOF is 17 and 15 for Xgr and Xmg> respectively 



4 Results - Comparison with observations - 

In this section we demonstrate the constraint on the modified gravity parameters /x and e 
using the observational data. We first use the surface mass density Xs(r_|_) and its logarithmic 
slope dlnSs/dlnrj_, obtained by Umetsu et al. through the accurate strong- and weak- 
lensing measurements [76]. We define the chi-squared for the surface mass density by 



X 2 {n,e,M vh ,c y 



n 



4 heo (^) 



Gov" 1 



->obs 



(rjj)) 



(4.1) 



where Eg bs (rj_j) and Cov^ 1 are the observational surface mass density and the inverse matrix 
of the la covariance matrix, respectively, at the i-th projected radius rj_j, which are obtained 
by stacking individual full surface mass density of four high-mass clusters (A1689, A1703, 
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A370, and C10024+17) at an average redshift (*,) = 0.32 [75, 76], and Eg heo (r_Lj) is the 
theoretically predicted surface mass density. 

For comparison, we also compute the chi-squared for the logarithmic slope d In Sg/ci In r± , 
which can be defined in the similar way to (4.1), 

y?( llfr M ) ~ V [d^r°(r ±l )/dln ri _ - din E° s bs (r ±t )/rfln r ± ] 2 
X^,e,c™,Mvit)-2_, [Ad\nX s (r ±i )/dlnr ± ] 2 ' [ ' 

where d\nT l f co (r u )/d\nr ± and d\nY>f s {r ±i ) / dhi r± are the theoretical value and the ob- 
servational value, respectively, and Ad In Xs(r_u)/d hi r ± is the 1-sigma error of the data (see 
the right panel of figure 1). 

The left panel of figure 11 shows the contour of Ax 2 , to show the relative confidence 
level of ji and e with respect to their best-fit values 

A% (Z^) *0 = X (M) ^' -^vir,local mini Cvir.local min) X (/^minj ^mim -^vir,min; Cvir,min) . (4.3) 

In (4.3), /x mm , e m i n , M v - lv ^ m \ a , and Cvir,min denote their best-fit values which globally minimize 
the value of x 2 > while M v ; ri i oca i m j n and c v i ri i OC ai min are the values that locally minimize the 
X 2 for a given set of values of \i and e, where we consider the range 

10 13 M Q < Mvir < 1O 16 M , (4.4) 

0.01 < c vir < 40, (4.5) 

and 

0.1 < e < 10, (4.6) 

-5.0 < n < 5.0. (4.7) 

Our results does not depend on this choice significantly. In table 1 we listed the best-fit value 
of the chi-squared of (4.1) Xmg = X 2 (^mm, e m in, M vir)min , c virimin ), for each halo model. 

In the left panel of figure 11, the region sandwiched between the solid curves is the 1- 
sigma allowed region with the chi-squared of the surface mass density (4.1), while the region 
between the dashed curves is the same but with the chi-squared of the logarithmic slope 
(4.2). Here we assumed the NFW profile. The small cross in this figure marks the value of 
the original galileon model, \i = 0.26 and e = 0.53 at the redshift 0.32 in the left panel. 

Note that we adopted the NFW profile in figure 11, while the panels in figure 12 assumed 
the different profiles. The panels (a) and (b) adopted the gNFW profile with (7^,7;) = 
(0.7,2.8) and (1.1,3.1), respectively, while (c) and (d) adopted the Einasto profile with T = 
0.17 and 0.23, respectively. 

We also investigate the constraint from the data of the differential surface mass density 
AS + (ri), obtained by Oguri et al. [77]. We define the chi-squared for AS + (rj_) by 

[A£V heo (rxi)-A£? s (r ±i ))] 2 



4 <b S ( r ^ 



where AS^ s (rj_j) and <r b s (r_u) are the observed differential surface mass density and the la 
error for the i-th projected radius rj_j in [77], and AS^ eo (r^j) is the theoretical differential 
surface mass density. The right panel of figure 11 is the contour of A% 2 for the differential 
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Figure 11. (a) 1-sigma confidence contour A\ 2 = 2.3 on the /i-e plane with the surface mass density 
Ss(?'_l) (solid curve), and with the logarithmic slope of the surface mass density din E s //dln rj_ 
(dashed curve). The region sandwiched between a pair of two curves is statistically allowed at 1- 
sigma level. The point marked by the cross corresponds to the original galileon model (see appendix 
B). (b) 1-sigma confidence contour Ax 2 — 2.3 on [i-e plane with the differential surface mass density 
A£ + (rjJ. In these panels, the NFW profile is assumed. 

surface mass density, denned in the same way as (4.3). The region sandwiched between the 
solid curves is the 1-sigma allowed region, and the point marked by the cross is the original 
galileon model, ji = 0.19 and e = 0.43 at the redshift z = 0.47. In the right panel of figure 
11, we adopted the NFW profile, and the best-fit value of the chi squared is noted in table 2 

2 

as Xmg- 

The panels in figure 13 are the same as the right panel of figure 11, but with adopting 
the different profiles. The panels (a) and (b) of figure 13 assume the gNFW profile with (7,,, 
7/) = (0.7,2.8) and (1.1,3.1), while (c) and (d) adopt the Einasto profile with V = 0.17 and 

0. 23, respectively. Note that the panel (d) indicates that the modified gravity is favored than 
Newtonian gravity when the Einasto profile with large V is assumed. 

5 Discussion 

Our main result is symbolically expressed by figure 11. From figure 11, Newtonian gravity, 

1. e., /i = or e > 1, is favored, and the original galileon model is in the 1-sigma allowed 
region, which is not rejected. One can see that the data put a constraint |^| ^ 1 when e^, 1, 
but the constraint on fj, becomes weaker as e becomes larger. For instance, —0.4 ?S /i ^ 0.6 
for e = 0.1 and -1.5 < fi < 2.5 for e = 1.0. 

The Vainshtein radius ry depends on e and M v ; r , from eq. (2.19). The method in 
the present paper cannot constrain the model in which the Vainshtein radius ry is larger 
than the cluster scale because the effect of the modified gravity becomes very weak inside 
the Vainshtein radius. The data we used is limited by r £ 3 /i _1 Mpc for the surface mass 
density, and by r ^ 5 /i _1 Mpc for the differential surface mass density. These facts set a 
limit of constraining /i and e. On the other hand, when the Vainshtein radius ry is smaller 
than the cluster scales r v j r , the effect of the modified gravity appears inside the cluster scale. 
Therefore, we can put a rather tight constraint on fj, for small e. 

In tables 1 and 2, respectively, the best-fit value of the chi-squared is listed for each halo 
model, using the surface mass density and the differential surface mass density. Xmg ^ s ^ ne 
best-fit value allowing the modified gravity model. Note that Xmg ^ s improved, compared 
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Figure 12. Each panel is the same as the left panel of figure 11, but adopted the different halo profile. 
The panels (a) and (b) adopted the gNFW profile with (7 S ,7;) = (0.7,2.8) and (j s ,ji) = (1.1,3.1), 
while (c) and (d) adopted the Einasto profile with T = 0.17 and T = 0.23, respectively. 
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Figure 13. Same figure as figure 12 but with the differential surface mass density AY, + (r±). 
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with Xgr' which is the best-fit value within Newtonian gravity. We see slight improvement 
for the NFW profile and the gNFW profile. This reflects that Newtonian gravity is favored 
from the data as long as the halo profile is the NFW profile or the gNFW profile. 

Another point is that the best-fit value is significantly improved by allowing the modified 
gravity for the Einasto profile with T = 0.23 and T = 0.3. This can be explained as follows 
by using figures 14 and 15, which show how the fitting is improved by the modified gravity 
model. Figure 14 shows that the fitting of Xs(r^) is improved at large radii by the modified 
gravity model. Figure 15 shows that the fitting of AS_|_(rj_) is improved at large radii as 
well as at smaller radii. Thus the modified gravity model is favored if the halo follows the 
Einasto profile with large value of T. This is the reason why the modified gravity is favored 
in the panel (d) of figure 13. 

The above results rely mostly on the behavior of the theoretical curves in the outer 
region. Let us check the validity of our theoretical modeling in the outer region. The 
behavior in the outer region is determined by the halo density profile which we choose and 
the modified gravity effect of the scalar field. When the density profile is fixed, the modified 
gravity's effect of the scalar field is properly taken into account by using the expression 
(2.13). In the present paper, however, we assumed the halo density profiles irrespective of 
the modified gravity model. This point will need to be tested carefully, though there exists a 
supporting evidence in reference [46, 51]. In the modeling of the halo density profile, we only 
considered the contribution of the cluster halo profile itself (1-halo term), and we neglect the 
neighboring halos (2-halo term), which will be only important at large radius [94]. 




Figure 14. This figure compares the observed surface mass density Ss( r ±) an d the theoretical curve 
assuming the Einasto profile with T — 0.23 (left panel) and T = 0.3 (right panel), respectively. In 
each panel, the blue dot and short-dashed curve is the best-fit curve within the Newtonian gravity, 
while the green dashed curve is the best-fit curve of the modified gravity. The solid curve assumes 
the NFW profile within the Newtonian gravity. The left (right) panel adopted T = 0.23(0.3), and the 
best-fit parameters for the modified gravity model are /i = 5.0 and e = 4.5 (fi = 3.4 and e = 0.69), 
respectively. 



6 Summary and Conclusion 

We investigated a constraint on the subclass of the generalized galileon model endowed with 
the Vainshtein mechanism, which includes a wide class of scalar-tensor gravity theories. 
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Figure 15. Same figure as figure 14, but for the differential surface mass density A£ + (rj_). The 
best-fit parameter for the modified gravity is fi = 5.0 and e = 1.8 (/i = 5.0 and e = 1.1) in the left 
(right) panel, respectively. 



The generalized galileon model we considered is characterized by the two parameters, \i 
and e, which describe the amplitude and the typical scale of the modification of gravity 
of the Vainshtein mechanism, respectively. We have investigated how the transition in the 
Vainshtein mechanism in the generalized galileon model appears due to a scalar degree of 
freedom, which leads to a modified relation between the lensing potential and the density 
perturbations of matter. This effect provides us with a unique opportunity for testing gravity, 
using observational quantities measured through the gravitational lensing. 

In the present paper, we focused our investigation on the density profile of cluster halo 
observed with gravitational lensing. Utilizing the observational data of the surface mass 
density in [76] and the differential surface mass density in [77], we obtained the constraint 
on the parameters \x and e for the first time, assuming the NFW profile, the gNFW profile 
and the Einasto profile as the halo density profile. We obtained the constraint, for instance, 
—0.4 ;$ fi £ 0.6 for e = 0.1, and — 1.5 ^ fi ^ 2.5 for e = 1.0, respectively, at 1-sigma confidence 
level, assuming the NFW profile (see figure 11). The constraint on \x is better for small e, 
but the constraint on /x becomes rather weak as e becomes larger. Newtonian gravity is 
favored when the NFW profile is assumed. This result is not altered when the gNFW profile 
is assumed (see section 4). However, the modified gravity is favored than Newtonian gravity 
when the Einasto profile with large T is assumed (see section 5). The constraint obtained with 
our method is not very tight. The original galileon model is not excluded in our test. But this 
method provides us with a unique chance to test the gravity theory on halo scales with cluster 
surveys, such as Subaru Hyper Suprime-Cam survey [98], CLASH [100], LoCuSS [101], and 
XXM-XXL [102]. 

In the present paper, we assumed the well-known density profiles of halo, which should 
be carefully tested with N-body simulations. In reference [46], properties of halos in the 
DGP modified gravity model are investigated. The DGP model possesses the Vainshtein 
mechanism similar to the galileon model. 3 From the N-body simulations in reference [46], it 
is suggested that the NFW profile describes the density profile of halos in the DGP model 

3 In the DGP model the relation connecting the lensing potential and the density perturbation is not 
altered, therefore it is difficult to obtain a constraint with the method in the present paper. 
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though the parameters of the profiles are scaled. This supports our assumption. However, 
in order to make the constraint robust, we need to check the modeling of the halo profile 
in future. As we considered the simple halo model in our analysis, our modeling of halo is 
valid only within the radius smaller than a few Mpc. In order to use the observational data 
including larger radius [94], we need a more sophisticated formulation, which is also left as a 
future problem. 
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A The coefficients in the perturbation equations 

In this appendix, we present the general expressions of the coefficients in the perturbation 
equations in a subclass of general second-order scalar-tensor theory with the action (2.1). 
Following [36], we introduce 

a °-jp + h- f - 2 ^t-^ip> (A - 1} 



where 



A 



H ' 



'G X 



o 



2H ' 



£ = 2XK X -K + 6HX(j)Gx - 2XG<p - 3H 2 F - 3iT<^, (A.2) 
V = K- 2X{G <t> + 4>G X ) + (3H 2 + 2H)F + (4> + 2H^)F <t> + 2XF H , 
B = -X4>G X + HF + ^/2, 



and F^ = dF((j))/d(j), Gx = dG/dX and G^ = dG/dcj). Finally, we have the following 
expressions for the coefficients in section 2, 

F<t> , c 2XG X -F 4> 2(A 1 +A 2 )H 2 B Q H<j> 

OL = -^-<P: C = — 0, C = = 5 A = — , 

F 2F [3Xcp 

F ^ + (A +A \ M 



P = -(a g+ aM + ^ 1 + A 2 )^Y4-. (A.3) 
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B Galileon model 



We consider the galileon model in curved spacetime with the action [29], 



S = d^xyf—g 



M 2 r 2 
2 Mpi r 



(B.l) 



which is obtained by taking F(4>) = M|j, K(X) = -X and G(X) = (r 2 /M Pl )X, where r c 
is the parameter. This model admits a late-time de-Sitter attractor in a flat FRW universe. 
The modified Friedmann equation on the attractor can be written as 



H(a)y 1 



Hn / 2 



n a~ 3 + J (n a~ 3 ) 2 + 4(1 - n c 



(B.2) 



The parameter r c is related to the cosmological parameters by r c = 1/(54(1 — Qq)) 1 ^}!^ 1 . 
In this galileon model, a, £, X 2 and (3 in (A. 3) are given by the background expansion 
history, as follows, 



_ 2XGx . 2 _ 4>Gx 



a = 0, t = M-{X4>G x , C = y 

P<P P 

l3 = -l + 2G x {ip + 2H^)-2M- 2 X 2 G 2 x . (B.3) 

For this galileon model along the attractor solution, it is useful to rewrite the combinations 
and X 2 ( in terms of the matter density parameter Q m = p m (a)/3Mp 1 i/ 2 (a) [33, 54], 

(1 — fi m )(2 — £l m ) 2 A ( 2 — J7 m \ , . 



In this model, we have /i = 0.26 and e = 0.53 (jt = 0.19 and e = 0.43) at z = 0.32 (z = 0.47), 
which corresponds to the mean redshift of the clusters to measure the surface mass density 
by Umetsu et al. in [75, 76] (the differential surface mass density by Oguri et al. in [77]), 
respectively. 
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